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Abstract 

We examine an Unruh-DeWitt particle detector coupled to a scalar field in 
three-dimensional curved spacetime. We first obtain a regulator-free expression 
for the transition probability in an arbitrary Hadamard state, working within 
first-order perturbation theory and assuming smooth switching, and we show that 
both the transition probability and the instantaneous transition rate remain well 
defined in the sharp switching limit. We then analyse a detector coupled to a 
massless conformally coupled field in the Hartle-Hawking vacua on the Bahados- 
Teitelboim-Zanelli black hole, under both transparent and reflective boundary 
conditions at the infinity. A selection of stationary and freely-falling detector 
trajectories are examined, including the co-rotating trajectories, for which the 
response is shown to be thermal. Analytic results in a number of asymptotic 
regimes, including those of large and small mass, are complemented by numerical 
results in the interpolating regimes. The boundary condition at infinity is seen to 
have a significant effect on the transition rate. 
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1 Introduction 



The Unruh-DeWitt model for a particle detector [TJ |2] is an important tool for probing 
the physics of quantum fields wherever noninertial observers or curved backgrounds are 
present. In such cases there is often no distinguished notion of a "particle," analogous to 
the plane-wave modes in Minkowski space, but an operational meaning can be attached 
to the particle concept by analysing the transitions between the energy levels of a detec- 
tor coupled to the field: upwards (respectively downwards) transitions can be interpreted 
as due to absorption (emission) of field quanta, or particles. The best-known applica- 
tions of this procedure are those for which the spectrum of transitions is thermal, which 
is the case for uniformly linearly accelerated detectors in Minkowski space [H [21 EJ H] , 
detectors at rest in the exterior Schwarzschild black hole spacetime El E] , and inertial 
detectors in de Sitter space [5|. 

With the Unruh-DeWitt detector, the fundamental quantity of interest is the prob- 
ability of a transition between the energy eigenstates. In the framework of first order 
perturbation theory this probability is proportional to the response function, given by a 
Fourier transform of the Wightman distribution of the quantum field over the detector's 
worldline, weighted by a switching function that specifies how the interaction is turned 
on and off [91 [10] . The response function is mathematically well defined provided the 
state of the quantum field is regular in the Hadamard sense [TT1 fT2] and the detector 
is switched on and off smoothly [131 HH US HE] • Physically, the response function then 
gives the probability for the detector to have completed a quantum jump by the time 
the interaction with the field has ceased. 

A related quantity of interest is the transition rate, which can be defined as the 
derivative of the transition probability with respect to the total detection time and is 
observationally meaningful in terms of consequent measurements in identical ensembles 
of detectors [T7]. There are technical subtleties in isolating in the transition rate the 
effects that are merely due to the switch-on and switch-off and the effects that are 
genuinely due to the acceleration and to the quantum state of the field [HI [191 EH 
EH [22l E3]; however, a satisfactory treatment within first order perturbation theory 
is to start with a smoothly-switched detector and take a controlled sharp switching 
limit pm I2"lf 125]. In particular, in three-dimensional Minkowski space, with a massless 
scalar field in the Minkowski vacuum, this procedure yields a finite result both for the 
transition probability and the transition rate even in the sharp switching limit [25] . 

In this paper we consider a detector in three-dimensional curved spacetime. In 
the first part of the paper we investigate a detector that is coupled to a scalar field 
in an arbitrary Hadamard state in an arbitrary spacetime. We give for the transition 
probability an expression that involves no distributional integrals, and we show that both 
the transition probability and the transition rate remain finite in the sharp switching 
limit. In the special case of a massless scalar field in the Minkowski vacuum of Minkowski 
spacetime, we recover the results of [25] . 

In the second part of the paper we apply the transition rate formula to a detec- 
tor in the (2 + 1) -dimensional Banados-Teitelboim-Zanelli (BTZ) black hole space- 
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time [261 127], f° r a massless conformally coupled scalar field in a Hartle-Hawking vacuum 
state [HI El ESI EH] with transparent, Dirichlet or Neumann boundary conditions at the 
asymptotically anti-de Sitter infinity. We first consider a stationary detector outside the 
hole, switched on in the asymptotic past. When the detector is co- rotating with the black 
hole horizon, we verify that the transition rate is thermal in the co-rotating local Hawk- 
ing temperature, in the sense of the Kubo-Martin-Schwinger (KMS) property [30], [3T] , as 
is to be expected from general properties of the Hartle-Hawking state P EJ ESI E2] and 
also from global embedding space (GEMS) considerations [331 EH [33 EE|- We obtain 
analytic results in a number of asymptotic regimes, including those of large and small 
mass, and we provide numerical results in the interpolating regimes. A static detector 
outside a nonrotating black hole is included as a special case. When the detector is 
stationary but its angular velocity differs from that of the horizon, we find that the 
transition rate breaks the KMS property already to quadratic order in the difference of 
the angular velocities of the horizon and the detector. 

We also consider a detector that falls into a nonrotating BTZ hole along a radial 
geodesic. We analyse the time evolution of the transition rate, obtaining analytic results 
for large black hole mass and large detector energy gap and numerical results in other 
regimes. We find no evidence of thermality in the transition rate, not even near the 
moment of maximum distance from the black hole, and we trace this phenomenon to 
the non-thermal response of an inertial detector in the corresponding vacuum state in 
pure anti-de Sitter space, as predicted from GEMS considerations [331 EH ESI 136] . 

We find that the boundary condition at the infinity has a significant effect on the 
transition rate for all of our detector trajectories. Typically, the Dirichlet condition 
gives a transition rate that varies least rapidly as a function of the detector's energy 
gap and the total detection time, owing to partial cancellations between terms in the 
Dirichlet Wightman function. Similar cancellations do not occur for the Neumann or 
transparent boundary condition. 

We begin in Section El with a brief review of the Unruh-DeWitt detector model, re- 
calling in particular how the distributional character of the Wightman function needs to 
be addressed prior to taking a sharp switching limit. In Section El we consider an arbi- 
trary Hadamard state in an arbitrary three-dimensional spacetime, writing the transition 
probability without distributional integrals and obtaining the formulas for the transi- 
tion probability and the transition rate in the sharp switching limit. Section H] briefly 
reviews relevant properties of the the BTZ black hole and the Hartle-Hawking vacua for 
a massless conformally coupled scalar field. A stationary detector is considered in Sec- 
tion [5] and a freely-falling detector in Section |6j Section [7] presents concluding remarks. 
Technical steps in contour integral analyses and asymptotic expansions are delegated to 
six appendices. 

Our metric signature is ( — h+), and we use units in which c — h — 1. Spacetime 
points are denoted by sans-serif letters. 0(x) denotes a quantity for which 0(x)/x is 
bounded as x — > 0, o(x) denotes a quantity for which o(x)/x — > as x — > 0, 0(1) 
denotes a quantity that is bounded in the limit under consideration, and o(l) denotes a 
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quantity that vanishes in the limit under consideration. 

2 Unruh-DeWitt detector 

In this section we briefly review the Unruh-DeWitt detector coupled to a scalar field, 
treated in first-order perturbation theory in the coupling (TJ |2] . 

The detector is an idealised atom with two energy levels, denoted by \0)d and \E)a, 
with the respective energy eigenvalues and E, where E may be positive or negative. 
The detector is spatially pointlike and moves on a timelike worldline x(r), parametrised 
by the detector's proper time r, in a spacetime of dimension two or higher. The detector 
interacts with a free real scalar field 0, of arbitrary mass and curvature coupling, by the 
interaction Hamiltonian 

Hi* = c X (t)^t)<P{x(t)) , (2.1) 

where c is a small coupling constant and /j,(t) is the atom's monopole moment operator. 
X is the switching function, positive during the interaction and vanishing elsewhere; this 
function specifies how the detector is switched on and off. We assume x to be smooth 
and of compact support, and we assume the trajectory to be smooth on the support 
of X- 

Before the interaction begins, we assume the detector to be in the state |0)d and the 
field to be in a state \^), and we assume |W) to be regular in the sense that it satisfies the 
Hadamard property [12]. The detector- field system is hence initially in the composite 
state \0)d <E> \^f). We are interested in the probability for the detector to be found in 
the state \E)d after the interaction has ceased, regardless of the final state of the field. 
Working in first order perturbation theory, this probability factorises as [51 dD] 

P(E) = c 2 \ d (0\fi(0)\E) d \ 2 F(E) , (2.2) 

where the response function F(E) encodes the information about the detector's trajec- 
tory, the initial state of the field and the way the detector has been switched on and off. 
J-{E) can be expressed as [T9] 

/oo /*oo 
du X (u) / dsx(u- s)e- iEs W e (u,u- s) , (2.3) 
oo Jo 

where W e (u, u — s) is a one-parameter family of functions that converge to the pull-back 
of the Wightman distribution to the detector's worldline [HI Q2], [EH dZ] . The factors 
in front of J~(E) in (12. 2p depend only on the internal structure of the detector and we 
shall from now on drop them, referring to F{E) as the transition probability. 

In summary, the response function F{E) answers the question 'What is the proba- 
bility of the detector to be observed in the state \E) d after the interaction has ceased?". 
Given F{E), we may define the detector's transition rate as the derivative of ^(E) with 
respect to the total detection time [HI [23] . Both JF{E) and the transition rate depend 
not only on the initial state of the field and the detector's trajectory but also on the 
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switching, and their behaviour in the sharp switching limit depends sensitively on the 
spacetime dimension [25] . 

Working directly with the expression ( 12 .3p for F(E) is cumbersome because the limit 
e — > 0+ may not necessarily be taken pointwise under the integral [171 HH EHl EH [22J 
[231 EU [25] . I n Section [3] we shall address this issue in three spacetime dimensions. 

3 Transition probability and transition rate in three 
spacetime dimensions 

In this section we specialise to three spacetime dimensions. We first rewrite the response 
function ( 12. 3 p in a form in which the regulator e does not appear. We then take the 
sharp switching limit and show that both the transition probability and the transition 
rate remain well defined in this limit. We follow closely the procedure developed in [TTj 
[231 E2 E5]. 

3.1 Hadamard form of W e 

In a three-dimensional spacetime, the Wightman distribution W(x,x') of a real scalar 
field in a Hadamard state can be represented by a family of functions with the short 
distance form [T2] 



where e is a positive parameter, <t(x, x') is the squared geodesic distance between x and x', 
cr e (x, x') := a(x, x') +2ie [T(x) — T(x')] + e 2 and T is any globally-defined future-increasing 
C°° function. The branch of the square root is such that the e — > 0+ limit of the square 
root is positive when <t(x, x') > [HI [T2]. Here U(x,x') and H(x,x') are symmetric 
biscalars that possess expansions of the form 



W £ (x,x') 



1 t/(x,x') 



+ 



H(x,x') 
V2 



(3.1) 




oo 




(3.2a) 



n=0 



oo 




(3.2b) 



n=0 



where the coefficients U n (x,x') satisfy the recursion relations 



(n + l)(2n + l)U n+1 + (2n + l)f/ n+ i ;M a^ - ( 
+ - m 2 - £R) U n = 0, n = 0,1,2,..., 



(2n + l)f/ n+1 A- 1 / 2 A 1 / 2 ;M ^ 



(3.3) 



with the boundary condition 



U = A 1 / 2 , 



(3.4) 
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and the coefficients H n (x,x') satisfy the recursion relations 

(n + l)(2n + 3)# n+1 + 2(n + l)H n+1 ,^ - 2(n + LU^A^A 1 / 2 ;^ 

+ {u x -m 2 -£R)H n = 0, n = 0,1,2,..., (3.5) 

where a = |cr, A(x, x') is the Van Vleck determinant, m is the mass and £ is the 
curvature coupling parameter [T2"] . 

The ze prescription in fl3.ll) defines the singular part of W(x, x'): the action of the 
Wightman distribution is obtained by integrating W e (x, x') against test functions and 
taking the limit e — > 0+ as in (12.31) . This limit can be shown to be independent of the 
choice of global time function T [131 H2 H3 EE] . 



3.2 Transition probability without it regulator 

To evaluate the e — > 0+ limit in (|2.3j) . the main issue is at s = 0, where the Hadamard 
expansion (13. ip shows that the integrand develops a nonintegrable singularity as e — > 0+. 
We shall work under the assumption that any other singularities that the integrand 
develops as e — > + are integrable. This will be the case in our applications in Sections E] 
and[6j We note in passing that similar integrable singularities can occur in any spacetime 
dimension, and the four- dimensional results in [T7j should hence be understood to involve 
a similar assumption. 

We first split the s-integral in (12.31) into the subintervals (0, rj) and (77, 00) where 
77 = yfe. We then find the small e expansions of each integral and finally combine the 
results. 

Let i> denote the s G (77, 00) portion of the s-integral in (12. 3p . including the taking 
of the real part. Let Wq denote the pointwise limit of W e as e — > 0+. Replacing W e by 
Wq creates in 7> the error 

POO 

2 Re / ds x{u - s) e~ iEs [W e {u, u-s)- W (u, u - s)] 
1 f°° 

= — Re / ds X (u - s) e~ tEs U(u, u-s) 
2vr J v 

where a is evaluated at the pair (x, x') = (x(n),x(n-s)), AT := T(x(u)) -T(x(n-s)) 
and we recall that a < as the trajectory is timelike. The square brackets in (13. 6p have 
the same form as in the Minkowski analysis in [23] and obey similar estimates. The 
function U has the small s expansion 

U(u,u-s) = l + 0(s 2 ) (3.7) 

by virtue of fl3.2a|) . (13. 4p and the expansions a = —s 2 + 0(s 4 ) and A = l + 0(s 2 ). These 
observations show that the contribution from (13. 6p to 7> is 0(rj). We hence have 

poo 

I> = 2 ds X (u- s)Re[e- iEs W (u,u- s)] + O (rj) , (3.8) 
Jo 



+ 2ieAT + e 2 
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where we have extended the lower limit to at the expense of an error that is contained 
in the O (rj) term since (13.1 1) shows that taking the real part under the integral makes 
the integrand regular at s = 0. 

Let then /< denote the s G [0,7]) portion of the s-integral in ( I2.3p . including the 
taking of the real part. We have 



/< = 2 Re p ds X (u - s) e~ lEs W e (u, u - s) 
Jo 

U(u,u — s) 



— Re / ds x( u ~ 

Z7T 



_y/a + 2ieAT + e 2 



+ 



H(u, u — s) 



V2 



(3.9) 



Since H(u, u — s) is a regular function, its contribution to J< is O {rj). The contribution 
of the term involving U can be found using the estimates of [25] and the expansion ( 13. 7p . 
We find 

I < = x(u)/4 + 0( V ). (3.10) 

Combining ( 13. 8 p and (13 . 1 Oj) . and noting that their error estimates hold uniformly 
in u, it is immediate to take the e — > 0+ limit in (12. 3p . We find 

i poo poo poo 

F{E) = - du [ X {u)f +2 du X (u) / ds X (u - s) Re [e~ iEs W (u, u-s)]. 

4 J-oo J~oc JO 

(3.11) 

Note that the integrals in (13. lip are regular, at s = by the Hadamard short- distance 
behaviour of Wq, and at s > by our assumptions about the the singularity structure 
of Wo at timelike-separated points. 



3.3 Sharp switching limit and the transition rate 

Up to now we have assumed the switching function x to be smooth. When x approaches 
the characteristic function of the interval [tq, tq+t] in a sufficiently controlled fashion [T71 
EU [25], the integrands in (13. lip remain regular, and taking the sharp switching limit 
under the integral can be justified by dominated convergence. The transition probability 
takes the form 

F T (E) = — ^ + 2 f du f ° dsRe [e~ iEs W (u,u- s)] , (3.12) 
4 Jt Jo 

where At := t — t and the subscript r is included as a reminder of the dependence on 
the switch-off moment. Differentiation with respect to r shows that the transition rate 
is given by 

T T (E) = \ + 2 f ds Re [e~ iEs W {r, r - s)] . (3.13) 
4 Jo 
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Note that both A3. 12[) and (13.131) are well defined under our assumptions, and in the 
special case of a massless scalar field in the Minkowski vacuum they reduce to what was 
found in [23]. Spacetime curvature has hence not introduced new singularities in the 
sharp switching limit. 

Note also that the pre-integral term ^ in ( 13 . 1 3[) would have been missed if the limit 
e — > + had been taken naively under the integral in (12. 3p . Yet this term is essential: it 
was observed in [25] that without this term one would not recover the standard thermal 
response for a uniformly linearly accelerated detector in Minkowski vacuum [31 [37] , and 
we shall see in Section [5] that without this term we would not recover thermality for a 
co-rotating detector in the BTZ spacetime. 



4 Detector in the BTZ spacetime 

We now turn to a detector in the BTZ black hole spacetime [261 EZ], specialising to a 
massless conformally coupled scalar field in the Hartle-Hawking vacuum with transpar- 
ent or reflective boundary conditions. In this section we briefly recall relevant properties 
of the spacetime and the Wightman function. More detail can be found in the review 
in [21]. 

Recall first that three-dimensional anti-de Sitter spacetime AdS3 may be defined as 
the submanifold 

- f = -Tl - Tl + X 2 + Xl (4.1) 

in M 2,2 with coordinates (T 1; T 2 , X±, X 2 ) and metric 

dS 2 = -dT 2 - dT 2 + dX\ + dX 2 , (4.2) 

where I is a positive parameter of dimension length. The BTZ black hole is obtained as a 
quotient of an open region in AdS3 under a discrete isometry group ~ Z. Specialising to 
a nonextremal black hole, a set of coordinates that are adapted to the relevant isometries 
and cover the exterior region of the black hole are the BTZ coordinates (t, r, 0), defined 
in AdS3 by 



Xi = i\fa sinh y~j<i> - ~^tj , X 2 = ly/a-1 cosh y—t j 

T 1 = Va^cosh(^>-^) , T 2 = iV^Tsmh(^t- r —<)?), (4.3) 



where 



(r 2 -r 2 \ 

and the parameters r± satisfy |r_| < r + . The coordinate ranges covering the black hole 
exterior are r + < r < oo, — oo < t < oo and — oo < < oo, and the Z quotient is 
realised as the identification (£, r, 0) ~ (t,r,(j) + 2ir). The outer horizon is at r — > r + , 
and the asymptotically AdS3 infinity is at r — > oo. The metric takes the form 

ds 2 = -{N^fdt 2 + f~ 2 dr 2 + r 2 (d(f) + N+dt) 2 (4.5) 
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with 



f 



r 2 J 2 \ 1/2 J 
£2 4 r 2 / 2r 2 



where the mass M and the angular momentum J are given by 



M 



J = 2r i r_ 



(4.6) 



(4.7) 



and they satisfy | J| < M£. 

In a quantum state invariant under d^, the Wightman function on the black hole 
spacetime can be expressed as an image sum of the corresponding AdS3 Wightman 
function. If GU(x, x') denotes the AdS3 Wightman function, the BTZ Wightman function 
reads |2S] 

Gbtz(x,x / ) = ^G a (x,AV) (4.8) 

where Ax' denotes the action on x' of the group element (t, r, 0) i— > (t, r, + 2tt), and the 
notation suppresses the distinction between points on AdS3 and points on the quotient 
spacetime. The scalar field is assumed untwisted so that no additional phase factors 
appear in (14. 8p . 

We consider a massless, conformally coupled field, and the family of AdS3 Wightman 
functions [29] 



^(x,xO = - 



c 



'AX 2 fx,x' 



AX 2 (x,x') +4£ 2 



(4.9) 



where the parameter ( e {0, 1, —1} specifies whether the boundary condition at infin- 
ity is respectively transparent, Dirichlet or Neumann. Here AX 2 (x, x') is the squared 
geodesic distance between x and x' in the flat embedding spacetime IR 2,2 , given by 



AX 2 (x,x') := -(T\ - T[f - (T 2 - T' 2 ) 2 + {X 1 - X[f + (X 2 - X' 2 f 

and we have momentarily suppressed the ie prescription in (14.91) . 
With (l4~8j) and (l4~9j) . the transition rate (I3.13p takes the form 



(4.10) 



oo 

4 2vrV2 f-^ 



At/1 



ds Re 



-iEls 



(4.11) 



'AX 2 V AX '+ 2 ^ 

where we have introduced the dimensionless integration variable s := sji and written 
AX 2 :=AX 2 (x(r),A' 1 x(r-S))/(2f) 

= -1 + y/a(r)a(r') cosh[(r+/£) (0 - <\> - 27m) - [r./i 2 ) (t - t')] 



- ^(a(r) - 1) (a(r') - l) cosh[(r+/f) (t - t') - (r-/£) (</> - <p' - 2nn)] 



(4.12) 
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where the unprimed coordinates are evaluated at x(r) and the primed coordinates at 
x(t-£s). 

What remains is to specify the branches of the square roots in (14. lip . As s extends 
to a global time function in the relevant part of AdS3, the prescription (13.11) implies that 
the square roots in (14. lip are positive when the arguments are positive, and the square 
roots are analytically continued to negative values of the arguments by giving s a small 
negative imaginary part. 



5 Co-rotating detector in BTZ 

In this section we investigate the transition rate of a detector that is in the exterior 
region of the BTZ black hole and co-rotating with the horizon. As the detector is 
stationary, we take the switch-on to be in the asymptotic past. When the black hole is 
spinless, the detector is static. 



5.1 Transition rate and the KMS property 

The angular velocity of the horizon is given by (26J EH EH] 

n H = r_/(r + £), (5.1) 

and it has an operational meaning as the value that d0/dt takes on any timelike worldline 
that crosses the horizon. The worldline of a detector that is in the exterior region and 
rigidly co-rotating with the horizon reads 

£r + r r_r 
r = constant , t = = , <p = = , (5.2) 

\J r 2 — r\^/ r\ — r _ y r 2 — y r ,2 _ — r _ 

where the value of r specifies the radial location and r is the proper time. We have set 
the additive constants in t and <p to zero without loss of generality. 

Substituting (15. 2 p into (I4.12p and taking the switch-on to be in the asymptotic past, 
the transition rate (14. lip takes the form 



oo 

HE) = ] + — A E 



oo 

ds Re 



e ~iE£s 




1 



Q n — sinh 2 (Ss + mxr^jt) 

where 



sinh 2 (Hs + nixr^jC) 

(5.3) 



K n := (1 - a' 1 )' sinh 2 (wrr + /£), (5.4a) 
Q n :=K n + {a-l)-\ (5.4b) 
S := (2y/a- (5.4c) 
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a is given by (I4.4p . and we have dropped the subscript r from J 7 as the situation is 
stationary and the transition rate is independent of r. The square roots in (I5.3P are 
positive for positive values of the argument, and they are analytically continued to 
negative values of the argument by giving s a small negative imaginary part. Note that 
the integrand in (15. 3p has singularities at s > 0, at places where the quantity under a 
square root changes sign, but all of these singularities are integrable. 
We show in Appendix [A~J that (15. 3p can be written as 

e -/3E£/2 00 



F(E) = ^ cos(n/3£r_) 



x 



/ X 

n=—oo 
oo 







dy cos(^A) ( 1 - C ] 

\J K n + cosh y \JQn + cosh y 



or alternatively as 



HE) = n ,o FJ , r , ^— / 



2(e^ + l) 2tt J v/Qo + cosh^ 



-PEi/2 ~ r°° ii t \ 

+ Yl cos ( n PEr„) / dy cos(yPE£/7r) = S , 

71 n=1 Jo \ a/ K n + cosh 2 y V Qn + cosh 2 ?/ / 



(5.6) 



where 



P := 2ttV« - 1 • (5.7) 



It is evident from (15.51) or (15.61) that J 7 depends on E only via the dimensionless com- 
bination £/3E. It is further evident that J 7 has the KMS property [301 EI] 

P{E) = e~^ E P{-E). (5.8) 

The transition rate is hence thermal in the temperature (£/3) _1 . 

It can be verified that (£(3) _1 = (—goo)" 1 ^ 2 ^, where T = K /(27r), k is the surface 
gravity of the black hole with respect to the horizon-generating Killing vector d t + 
VLndcf,, and goo is the time-time component of the metric in coordinates adapted to 
the co-rotating observers. This means that the temperature (£[3)~ l of the detector 
response is the local Hawking temperature, obtained by renormalising the conventional 
Hawking temperature T by the Tolman redshift factor at the detector's location. This is 
the temperature one would have expected by general properties of the Hartle-Hawking 
state [61 [71 [28], including the periodicity of an appropriately-defined imaginary time 
coordinate [32], and also by GEMS considerations [331 EH ESI [36] . 

Note that the expressions (15.51) and (15. 6p contain both terms of ( 13. 1 3j) . as shown in 
Appendix [Al The pre-integral term j in (13.131) is hence essential for recovering ther- 

mality: in (15 .6p it can be regarded as having been grouped in the term ^(e fSE£ + 1) , 
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which gives the transition rate in pure ACIS3 with the transparent boundary condition. 
The superficial Fermi-Dirac appearance of this pure AdS3 term is a general feature of 
linearly-coupled scalar fields in odd spacetime dimensions [21 EBl E3 EE]- 



5.2 Asymptotic regimes 



We consider the behaviour of the transition rate (15. 6p in three asymptotic regimes. 

First, suppose r + — > 00 so that r_/r + and r/r + are fixed. Physically, this is the 
limit of a large black hole with fixed J/M, and the detector is assumed not to be close 
to the black hole horizon. Note that a and (3 remain fixed in this limit. It follows from 
( 14.41) and (15.41) that in (15.61) this is the limit in which K n and Q n with n > 1 are large. 
Assuming that E is fixed and nonzero, and using formula (lB.2al) in Appendix [B], we find 



HE) 



1 



Ce 



-/3EI/2 poo 



+ 



2{eP El + 1) 



2tt 



dy 



cos(y/3E£/ir) 
a/Qq + cosh 2 y 



x < Im 



\i/3El/(2ir) 



r 1 + 



i/3E£ 
2vr 



ipEi 
2vr 



+ 0(e 



(5.9) 



where the displayed next-to-leading term comes from the n = 1 term in (I5.6P and is 
of order e _7rr+//£ . The corresponding formula for E = can be obtained from formula 
(1B.2bl) in Appendix |B] and has a next-to-leading term of order r + e~ nr+ ^ e . 

Next, suppose that r + — > so that r„/r + and r/r + are again fixed. This is the 
limit of a small black hole. Note that a and are again fixed. The dominant behaviour 
comes now from the sum over n and can be estimated by the Riemann sum technique 
of Appendix O We find 



HE) 



g e -/3Ee/2 



+ JO 



00 />oo 

dv 



,vf3E£r\ fyPEC 
dy cos I cos I I x 



X 



asinh v 



V 7r 

2. 



cosh 2 y ) - C f — , _ v h cosh 2 y 



ofl] 



(5.10) 



The leading term is proportional to l/r + and it hence diverges in the limit of a small 
black hole. 

Finally, suppose that E — > ±00 with the other quantities fixed. The analysis of 
Appendix ITJ1 shows that each integral term in ( 15 .6p is oscillatory in E, with an envelope 
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that falls off as 1/y—E at E — > — oo but exponentially at E — > +00. Applying this 
estimate to the lowest few values of n in (15.61) should be a good estimate to the whole 
sum when r + /£ is large. We have not attempted to estimate the whole sum at E — > ±00 
when r + /£ is small. 



5.3 Numerical results 

We now turn to numerical evaluation of the transition rate ( 15. 6p . We are particularly in- 
terested in the interpolation between the asymptotic regimes identified in subsection 15. 21 

T (15. 6p depends on five independent variables. Two of these are the mass and the 
angular momentum of the black hole, encoded in the dimensionless parameters r + /£ and 
r_/£ The third is the location of the detector, entering T only in the dimensionless 
combination a (14.41) . The fourth is the detector's energy gap E, entering T only in the 
dimensionless combination (3E£ where (3 was given in (15. 7p . The last one is the discrete 
parameter ( G {0, 1, —1} which specifies the boundary condition at infinity. 

We plot J as a function of £f3E, grouping the plots in triplets where ( runs over 
its three values and the other three parameters are fixed. We proceed from large r + /£ 
towards small r + /£. 

In the regime r + /£ > 3, numerics confirms that the n > 1 terms in (15.61) are small. 
JF therefore depends on r + /£ and r_/£ significantly only through (3, that is, through the 
local temperature. The detector's location enters T in part via (3 (15.71) . but also via Qo 
in (15.61) . and the latter affects only the the boundary conditions ( = 1 and £ = — 1, in 
opposite directions. Plots for r + /£ = 10 are shown in Figured) 

As r + /£ decreases, the n = 1 term in (15.61) starts to become appreciable near r + /£ « 
1. The dependence on r_/£ is then no longer exclusively through 0, and the effect is 
largest for ( = and £ = — 1 but smaller for £ = 1, owing to a partial cancellation 
between the two terms under the integral in (15.61) for ( = 1. Plots for r + /£ = 1 are 




-30 -20 
3EI 



-a(r) =4 a(r) = 100 

(a) C = 



-50 -40 -30 -20 



-q(r) =4 q(r) = 10o| 

(b)C = l 



-50 -40 -30 -20 -10 

PEl 

\^a{r) =4 a{r) = 100 

(c)C = -l 



Figure 1: t as a function of f3E£ for r + /£ = 10 and r_j£ = 0, with a = 4 (solid) and 
a = 100 (dotted). Numerical evaluation from (15 .6p with n < 3. 
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J3E£ 



K=0 K=0.99| 

(a) C = 



(b)C = l 



(c)C 



K=0.99| 



-1 



Figure 2: J 7 as a function of f3E£ for r + /£ = 1 and a = 2, with k = (solid) and 
k = 0.99 (dotted) where k := r_/r + . Numerical evaluation from (15. 6p with n < 3. 




-40 -30 -20 -10 



-K=0 K=0.99| 



-40 -30 -20 -10 

PEl 

| K = K=0.99 



40 -30 -20 -10 



_PJ| 

-K=0 K=0.99| 



(a) C = (b) C = 1 (c) C = -1 

Figure 3: As in Figure [2] but for a = 100. 



shown in Figures [2] and [3j 

Asr + /£ decreases below 1, the next-to- leading asymptotic formula (15. 9p starts to be- 
come inaccurate near r + /£ = 0.3, as shown in Figure HJ although the partial cancellation 
between the two terms under the integral in (I5.6P and the similar partial cancellation in 
(15.91) moderates the effect for ( = I. At r + /£ = 0.1, shown in Figure El T is sensitive 
to changes in both r_/£ and a. When a ^> 1, the Q = — 1 curves in Figure [5] have 
approximately the same profile as the ( = curves but at twice the magnitude: from 
(I5.4ap and f!5.4bp we see that this indicates the regime where the n > 1 terms in (15. 6|) 
give the dominant contribution to T. 

As r + /£ decreases further, we enter the validity regime of the asymptotic formula 
(I5.10p . as shown in Figure [6] for r + /£ = 0.01. Note that again the ( = — 1 curve has 
approximately the same profile as the ( = curve but at twice the magnitude, indicating 
that the dominant contribution comes from the n > 1 terms in (15. 6p . 
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•4,6 
0.4 
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- Numerics ' ' ' ' Asymptotics | 
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10 -50 -40 -30 -20 -10 



(b) C = 1 



(c) C = -1 



Figure 4: j as a function of (3E£ for r + /£ = 0.3 and r_ji = 0.299, with a = 2. 
Solid curve shows numerical evaluation from (15. 6p with n < 3. Dotted curve shows the 
asymptotic large r + /£ approximation (15.91) . 




50 -40 -30 -20 -10 



-(2,0) (2,0.099)-"- (100,0) 
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-(2,0) (2,0.099) (100,0) 

(100,0.099) 
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Figure 5: T as a function of /3E£ for r+jl — 0.1, with selected values of the pair (a, r_/i) 
as shown in the legend. Numerical evaluation from (15.61) with n < 35. 
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-50 -40 -30 -20 -10 -50 -40 -30 -20 -10 -50 -40 -30 -20 -10 

3EI [3EI 3EI 

I Numerics • • • ' Asymptotics | | Numerics ■ ■ ■ ■ Asymptotics | | Numerics • • • ' Asymptotics | 

(a) C = (b) C = 1 (c) C = -1 

Figure 6: J as a function of jiEl for r + /£ = 0.01 and r_ = 0, with a = 4. Solid curve 
shows numerical evaluation from (15. 6p with n < 300. Dotted curve shows the asymptotic 
small r + /£ approximation f)5.10p . Qualitatively similar graphs ensue for r_/r + = 0.99. 



6 Radially in-falling detector in spinless BTZ 

In this section we consider a detector on a radially in-falling geodesic in a spinless BTZ 
spacetime. 



6.1 Transition rate 

Recall from Section H] that for a spinless hole r_ = and r + = Mi > 0, and the horizon 
is at r = r + . To begin with, we assume that at least part of the trajectory is in the 
exterior region, r > r + . Working in the exterior BTZ coordinates (I4.3p . the radial 
timelike geodesies take the form 



Mq cos f, 

tanf 



t = (£/VM) arctanh 



(6.1; 



where q > 1, </>o denotes the constant value of 0, and f is an affine parameter such that 
the proper time equals f£. The additive constants in f and t have been chosen so that 
r reaches its maximum value £y/Mq at f = with t = 0. 
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Substituting (16. ip in (14. lip and (I4.12p . we find that the transition rate is given by 



Me) = i/4 



+ 



2ttV2 



E 



ds Re 



-iEs 



-c 



■sj— 1 + K n cos f cos(r — s) + sin r sin(r — s) 

Q—iEs 

y/l + -ftT n cos f cos(f — s) + sin f sin(f — s) 



(6.2) 



where 

AT„ := 1 + 2g 2 sinh 2 (nvrv / M). (6.3) 

The detector is switched off at proper time r and switched on at proper time To = r-Ar, 
and we have written f := t/£, Af := At/£ and E := E£. The square roots in (16. 2p 
are positive when the arguments are positive, and they are analytically continued to 
negative values of the arguments by giving s a small negative imaginary part. 

Although the above derivation of (16. 2 p proceeded using the exterior BTZ coordinates, 
the result (16. 2p holds by analytic continuation even if the geodesic enters the black or 
white hole regions. The ranges of the parameters are — tt/2 < f — Af < f < 7r/2, so that 
the detector is switched on after emerging from the white hole singularity and switched 
off before hitting the black hole singularity. 



6.2 The n = term and KMS 

We write (16. 2 p as 

j~ t = j; n=0 + , (6.4) 

where J r ^ =0 consists of the n = term and J 7 ™^ consists of the the sum J^n^o- ^ e 
consider first ^"™ =0 . 

jrn=o gj ves transition rate of a detector on a geodesic in pure AdS3. ^™ =0 does 
not depend on M or q, and it depends on the switch-on and switch-off moments only 
through Af, that is, through the total detection time. Using (16. 2p and (16. 3p . we find 

T ^ ' 4 4W 

where E := El. As < Af < n, 06.51) is well defined. 

Numerical examination shows that T™ =Q does not satisfy the KMS condition. This is 
compatible with the embedding space discussion of [331 EH ESI EH], according to which 
a stationary detector in AdS3 should respond thermally only when its scalar proper 
acceleration exceeds l/l. 



sin (Es) cos(ES) 
sin(s/2) cos(s/2) 



(6.5) 
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The asymptotic behaviour of J 7 ™ at large positive and negative energies for fixed 
Af can be found by the method of Appendix [El We find 

• n/ Q(-E) 1 f cos(EAf) sin(£Af)\ , , ~„ 

rr° (E) = A-J. + ^-^A _ c^^j + o(i/n (6.6) 

where is the Heaviside step-function. 



6.3 The n^O terms and large M asymptotics 

We now turn to J 7 ™^ , which contains the dependence of T T on M and q. 

We consider J 7 ™^ in the limit of large M. We introduce a positive constant c G 
(0, 7r/2), and we assume that the switch-on and switch-off moments are separated from 
the initial and final singularities at least by proper time ci. In terms of f and Af, this 
means that we assume 

-tt/2 + c < f < tt/2 - c , < Af < f + tt/2 - c . (6.7) 

As iT n = K_ n , we can replace the sum Ylm^o * n by 2^^ =1 . Given (16. 7p . the 
expression cosf cos(f — s) is bounded below by a positive constant. Using f l6.3|) . this 
implies that the quantities under the n^O square roots in (16. 2p are dominated at large 
M by the term that involves K n , and we may write 



1 °° i /-Af 

7rv2cosr " Jo 



cos 



(£§)d§ / 1 C 



Vcos(f - 5) V\/l + /-/^n ^1 + /+/AT n 

(6.8) 



where 



u ■= sinfsin(f " g) j: 1 . (6.9) 



COS T COS T 



The large M expansion of J 7 ^ is then obtained by a binomial expansion of the square 
roots in (16. 8p at K n — > oo and using (16. 3p . The expansion is uniform in f and Af within 
the range (16 .7p . and by (I6.3P it is also uniform in q. The first few terms are 



. , n/ s 1 f Af cos(Es)ds 



7^^ 



\J2 cos f Jo ^cos(t 



^ + 



c/+ - /- 



+ 0(e~ 5WM ) . (6.10) 

For ( 7^ 1, the dominant contribution comes from the term proportional to (1 — () and 
is of order e~ n ^ . 
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6.4 Numerical results 



At large M, the dominant contribution to T T comes from f™ =Q (jfi.sp , which depends 
only on and At/£. Plots are shown in Figures [7] and |HJ When \E£\ is large, the 
oscillatory dependence on At/ £ shown in the plots is in agreement with the asymptotic 
formula (16. 6p . 

When M decreases, the contribution to T T from J 7 ™^ becomes significant. For 
M = 0.1, the terms shown in (I6.10p are still a good fit to the numerics provided both 
the switch-on and the switch-off are in the exterior region. For smaller M, the number 
of terms that need to be included in J 7 ™^ increases rapidly. A set of plots is shown in 
Figures 191 and fTUl for M = 1CT 4 with q = 100, taking the detector to be switched on at 
the moment where r reaches its maximum and following the detector over a significant 
fraction of its fall towards the horizon. J 7 ™^ turns out to be still insignificant at large 
negative E£, but it starts to become significant at E£ > —5, and its effect then depends 




Figure 7: j"^ =0 O as a function of E£ and At/£ for ( = 0. 




0.2 0.4 0.6 0.8 1 1.2 1.4 0.2 0.4 0.6 0.8 1 1.2 1.4 0.2 0.4 0.6 0.8 1 1.2 1.4 

At/1 Ar/i Ar/i 

(a) El = -100 (b) Ei = -5 (c) Et = 20 



Figure 8: (16.51) as a function Ar/£ for selected values of E£, with ( = (dashed 
line), ( = 1 (thick line) and ( = —1 (dotted line). 
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Figure 9: T T O with M = 10" 4 , q = 100, r = and Et = -5. Solid curve shows 
numerical evaluation from ( 16. 2p with 200 terms and dashed curve shows the individual 
n = term ^" =0 (16. 5ft . The horizon-crossing occurs outside the plotted range, at 
At/£ = arccos(O.Ol) « 1.56. 
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Figure 10: As in Figure M but with Ei = 20. 



strongly on the boundary condition parameter £, being the smallest for £ = 1. 

For fixed M, following the detector close to the future singularity numerically would 
pose two complications. First, an increasingly large number of terms would need to be 
included in J 7 ™^ . Second, the evaluation of the individual terms to sufficient accuracy 
would need to handle numerically integration over an integrable singularity in s. This 
singularity arises because the the quantity under the first square root in ( 16 .2p can change 
sign within the integration interval. We have not pursued this numerical problem. 



7 Conclusions 

In this paper we have investigated the response of an Unruh-DeWitt particle detector 
in three-dimensional curved spacetime. We first obtained a regulator-free expression for 
the transition probability of a detector coupled to a scalar field in an arbitrary Hadamard 
state, working within first-order perturbation theory and assuming that the detector is 
switched on and off smoothly, and we showed that both the total transition probability 
and the instantaneous transition rate remain well defined when the switching becomes 
sharp. In the special case of a detector in Minkowski space, coupled to the Minkowski 
vacuum of a massless scalar field, these results reduce to those found in [25]. The 
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results confirm that the sharp switching limit is qualitatively different between three 
and four spacetime dimensions even when the spacetime is curved: in four dimensions, 
the sharp switching limit yields a well-defined transition rate but a divergent transition 
probability fiT\fM\. 

We then specialised to a detector in the BTZ black hole spacetime, coupled to a 
massless conformally coupled scalar field in a Hartle-Hawking vacuum state with trans- 
parent, Dirichlet or Neumann boundary conditions at the infinity. For a stationary 
detector that is outside the horizon and co-rotating with the hole, and switched on 
in the asymptotic past, we verified that the transition rate is thermal in the sense of 
the KMS property, in the local Hawking temperature that is determined by the mass 
and angular momentum of the hole and by the Tolman redshift factor at the detector's 
location. This is the temperature that was to be expected by general properties of 
the Hartle-Hawking state E2], and by GEMS considerations [331 Ei ESI EE]- 
A static detector outside a nonrotating black hole was included as a special case. We 
obtained analytic results for the transition rate in a number of asymptotic regimes of the 
parameter space, including those of large and small black hole mass, and we provided 
numerical results in the interpolating regimes. We have not pursued in detail the case of 
a stationary detector whose angular velocity differs from that of the hole, but we shall 
show in Appendix [F] that the parameter space has at least some regimes in which the 
response of such a detector does not have the KMS property. 

We also considered a detector that falls into a nonrotating BTZ hole along a radial 
geodesic. As the trajectory is not stationary, the transition rate is not constant along 
the trajectory, and in particular the switch-on cannot be pushed to the asymptotic 
past since the trajectory originates at the white hole singularity. We obtained analytic 
results for the transition rate when the black hole mass is large, and we evaluated the 
transition rate numerically for small values of the black hole mass provided the switch- 
on and switch-off take place in the exterior. We found no parameter ranges where the 
transition rate would be approximately thermal in the sense of the KMS property, not 
even near the moment of maximum radius on a trajectory, and we traced the reasons for 
this to the properties of AdS3 geodesies that have been previously analysed from GEMS 
considerations (331 EH ESI EE] • Our expression for the transition rate as a countable sum 
remains valid after the detector crosses the horizon, but the sum becomes then more 
difficult to estimate analytically and more labour-intensive to evaluate numerically, and 
we have not pursued a detailed investigation of this regime. 

It would be interesting to compare our BTZ transition rates to those of a detector on 
similar trajectories in four-dimensional Schwarzschild space-time, where the Wightman 
function needs to be evaluated fully numerically. Some differences can be expected to 
arise from the different asymptotic infinities of BTZ and Schwarzschild: for example, 
an inertial detector in Schwarzschild should respond to the Hartle-Hawking vacuum 
approximately thermally in the asymptotically flat region. We leave this question subject 
to future work. 
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A Appendix: Derivation of (15.51) and (15.61) 



In this appendix we verify the passage from (15. 3 p to (|5.5p and (15. 6p . 

A.l n = term 

Let 

poo -iax J 

J(a,P):=Re/ (A.l) 
Jo v P — sinh x 

where o G K, ? > 0, and the square root is positive for positive argument and positive 
imaginary for negative argument. We shall show that 

J(a,0) = - 7rtanh 2 (7ra/2) , (A.2a) 

I(a, P) = e-/ 2 r ™s(ay)dy for P > . (A.2b) 

Jo a/ P + cosh 2 ?/ 



Applying (I A . 2 1) and (I A . 3 1) to the n = term in (I5.3P yields the corresponding terms in 



(1531) and (CT 

Suppose first P = 0. For P = 0, (lA.ip reduces to I(a, 0) = — / °° ^j^r dx, which 



evaluates to (lA.2ap [39|. 

We note in passing the relation 



71 , „-™/2 f°° cos(ay) dy 



Z (a ,0) = -- + e-^y o -±fiJL, (A.3) 

which follows by evaluating the integral in (IA.3[) [39J and using ( IA.2ap . Comparison of 
(1A.2bp and ( 1A.3|) shows that I(a, P) is not continuous at P = 0. 
Suppose then P > 0. We rewrite ( lA.ip as the contour integral 



J(a,P):=Re / \ , (A.4) 



i \fp~— sinh 2 2 
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where the contour C\ goes from z = to z = oo along the positive real axis, with a dip 
in the lower half-plane near the branch point z = arcsinh \f~P . The square root denotes 
the branch that is positive for small positive z. 

We deform C\ into the union of C2 and C3, where C2 goes from z = to z = —in/ 2 
along the negative imaginary axis and C3 consists of the half-line z = y — in /2 with 
< y < 00. As the integrand has no singularities within the strip —it/2 < Imz < and 
falls off exponentially within this strip as Re z — > +00, the deformation does not change 
the value of the integral. The contribution from C2 is purely imaginary and vanishes on 
taking the real part. The contribution from C3 yields (IA.2bl) . 

A. 2 terms 

Let 

roc -iax J 

I b (a, P) := Re / = , (A.5) 

'P - sinh 2 (x + b) 







where aGK, P > 0, 6 G M and the square root is positive for positive argument and 
analytically continued to negative values of the argument by giving x a small negative 
imaginary part. 

We shall show that 

I b {a, P) + I- b (a, P) = 2 cos (ab) I(a, P) , (A.6) 



where I(a,P) is given in fTA~2b|) . Applying (IP)) with f lA~2bj) to the n ^ terms in 



yields the corresponding terms in (15.51) and (15.61) . 

For b = 0, (1A.6I) follows from (lA.2bj) . As both sides of ( 1A.6I) are even in b, it hence 



suffices to consider (1A.6|) for 6 > 0. 

Let 6 > 0. Changing the integration variable in f1A.5h to y = x + b yields 



I b (a, P) + I-b(a, P) = 2 cos (ab) Re 



dy 



yP — sinh 2 y 



Jo y P — sinh 2 y Jo y P — sinh 2 y 



where the branches of the square roots are as inherited from (1A.5|) : positive when the 



argument is positive and continued to negative argument by giving y a small negative 
imaginary part. Examination of the branches shows that the last two terms in (1A.7j) 
cancel on taking the real part, and using (lA.lj) in the first term leads to (]A.6|) . 



B Appendix: Derivation of (15.91) 



In this appendix we verify the asymptotic formula (15 .9|) . 
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Let 



J(a,P) := f 
Jo 



cos(ay) dy 
o \/~P + cosh 2 ?/ 



(B.l) 



where P > and a e R. Note from f[Ob|) that /(a, P) = e-™l 2 J(a, P) for P > 0. We 
shall show that as P — > oo with fixed a, J(a, P) has the asymptotic form 

J(a, P) = — ^= Im [(4P) ia/2 r(l + m/2)r(i - ia/2)] + O (P~ 3/2 ) for a ^ 0, (B.2a) 
av 7rP 



J(0,P) = -1= [ln(4P) +^(1) +0(P" 3 / 2 In P) 

Zy r 



(B.2b) 



where ^ is the digamma function [4"0] . 

Starting from (IB.lj) . writing cos(ay) = Re(e my ) and making the substitution y = hat, 
we find 



J(a,P) = 2Re J i 



t ia dt 



Vt A + BH 2 Jl 



t 4 + B 2 t 2 



OO „r 

2^6 p Re / 
p=0 ^ x 



t ia dt 



t 2P+i(t*+B 2 y +{1/2) 



(B.3) 



where B = \J\P + 2 and b p are the coefficients in the binomial expansion (1 + x) 1 ' 2 — 
Y^=o h P xP - As the P > terms in ( B) are C^P" 2 ^ 1 ) = C^P"^ 1 / 2 )) by dominated 
convergence, we have J(a,P) = J (a,P) + 0(P~ 3 / 2 ), where the substitution t = Bv in 
the p = terms gives 



9 / /»oo „,ia—l 

J J a) = - Re ( B ia 



B 



v ia 1 ^ v 



l/B a/1 + V< 



(B.4) 



When a 7^ 0, integrating (1B.4|) by parts and extending the lower limit of the integral 
to zero gives 



J (a,P) = — Im 
Ba 



B l 



v l+ia ^ v 



(1 + V 2 ) 



3/2 



0(P~ 2 ) 



(B.5) 



The integral in (1B.5|) may be evaluated by 

writing (1 + 1> 2 ) 3 ^ 2 = (r(3/2)) 1 J °° dn 1 ' 2 e -(i+*> 2 ):r anc [ interchanging the order of 
the integrals, with the result (lB.2aj) . When a = 0, similar manipulations lead to (1B.2bj) . 
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C Appendix: Derivation of (15.101) 

Let p > 0, q > 0, oGK and 76M. For n G Z, let K n := p 2 sinh 2 (ng), and define 

cos (7170) cos (ay) dy 



\/K n + cosh' 



(C.l) 



where we suppress the dependence of F n on p, q, a and 7. We shall show that the sum 
S := J2^=~oo Fn has the asymptotic form 

S = - f°dr [°° cos (nO cos(ay) dy + o(l) 
9 Jo Jo ^/ p 2 sinh 2 r + cosh 2 ?/ 9 

as g — >■ with the other parameters fixed. Note that the leading term in (10.2[) diverges 
as g — > 0. 
Let 

^, x / \ f°° cosjay) dy 

G(r) := cos( 7 r) / -===== == , (C.3) 

Jo v p z sinn r + cosh y 

where we suppress the dependence of G on a and 7. S then equals q~ l times the 
Riemann sum of G with the sampling points r = ng, n G Z. G is continuous, and from 
Appendix IB1 we see that |G(r)| is exponentially small as r — > ±00. The Riemann sum 
of G therefore converges to the integral of G as q — > 0. Noting finally that G is even, 
we recover (1C.2I) . 



D Appendix: Co-rotating response at E£ — > ±00 



In this appendix we analyse the individual terms in the co-rotating detector response 
(15.61) in the limit El — > ±00. These terms are of the form 



I( X ,a,P) :=cos( X a) e^ 2 J(a,P) 



(d.i; 



where x^^; a ^^ ; -P>0 and J(a,P) is given by (IB.ip . We shall show that when 
a — > ±00 with fixed x an d Pi P\Xi a ; P) has the asymptotic expansion 



'2y/ne na cos (%a) cos(aa — 7r/4) 
^/ a sinh(2a) 

2-\/7rcos (xa) cos(— aa — it/A) 



+ o(a 



yf— asinh(2a) 



-1/2 g-aTA 



+ o((-a)- 1/2 ) 



a +00, 



a — > —00, 



(D.2) 



where a = arcsinh \f~P . 
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Assuming a ^ and writing cos(ay) = Re(e*' a ' y ), we start by rewriting J(a,P) from 
flRTj) as 

J(a,P) = Re / - \ , (D.3) 

JCi \/P + cosh 1/ 



where the contour C\ consists of the positive imaginary axis travelled downwards and 
the positive real axis travelled rightwards. The contribution from the imaginary axis 
vanishes on taking the real part. 

Writing P = sinh 2 a where a > and factorising the quantity under the square root 
in flD.3p . we obtain 

J(a, P) = Re / - U . (D.4) 

</Ci \ sinhfa + y — in/2) sinhfa — y + m 12) 



Ci v sinh(a + y — in/2) sinh(a — y + m/2) 

The branch points of the integrand in (ID. 41) are at y = ±ct + in(n + A n £ Z. We may 
deform C\ into the contour C2 that comes down from a + zoo at Re y = a, passing the 
branch points from the left, encircles the branch point at y = a + m/2 counterclockwise, 
and finally goes back up to a + zoo at Re y = a but now passing the branch points from 
the right. Changing the integration variable by y = a + m /2 + iu, we then have 

J(a,P) = e-H-/2 Re (, e Ha| f e-^ du \ ^ ^ g) 

\ Jc 3 y—i sin(w) sinh(2a + iu) J 

where contour C3 comes from u = +00 to u — on the upper lip of the positive u 
axis, encircles u = counterclockwise and goes back to u — +00 on the lower lip of 
the positive u axis. The square root is positive at u = ir/2 on the upper lip and it is 
analytically continued to the rest of C3. 

We now note that sinh(2a + iu) = sinh(2a) cos(w) + i cosh(2a) sin(w), and that the 
modulus of this expression is bounded below by sinh(2a). In (1D.5j) . the contribution from 
the two intervals in which it/2 < u < it is therefore bounded above by e~^ w / \J sinh(2a) 
times a numerical constant, and the contribution from the two intervals in which inn < 
u < (n + l)vr, n = 1,2,..., is bounded above by e~\ a \ 7r l- n+ ( 1 / 2 ^ / sinh(2a) times a 
numerical constant. The sum of all of these contributions is hence 0(e~^ n ). In the 
remaining contribution, coming from the two intervals in which < u < tt/2, we combine 
the upper and lower lips and change the integration variable to w = \a\u. This gives 



2e -|a|7r/2 

J(a,P) = 

\/\a\ 



kk/2 e -w dw 



x Re ( e i(a|a| "" /4) . 

sin(u>/|a|) [sinh(2a) cos(to/|a|) + i cosh(2a) sin(w/|a|)] 
0(e- |a| ") , (D.6) 
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where the the square root denotes the branch that is positive in the limit w — > + . 

By Jordan's lemma, the modulus of the integrand in flD.6j) is bounded from above 
in the range of integration by the function g(w) := 2 sinh(2a j w ~ 1 ^ 2 e ~" '• As 9i w ) is 
integrable over < w < oo and independent of a, dominated convergence guarantees 
that when \a\ — > oo, the limit in the integrand in (ID.6[) can be taken under the integral. 
The integral that ensues in the limit is elementary, and we obtain 

J(a, P) = 2 ^ FeHa '; /2cOS(Q|a| - 7r/4) + o(\a\-^ e-H^) . (D.7) 
y\a\ sinh(2o;) 

flTX2]) then follows by substituting f lTX7j) in flDvLj) . 



E Appendix: Derivation of (16.61) 



In this appendix we verify the asymptotic expansions 

cos (f3x) sin(m/3) 



dx 



dx 



+ 0(/3" 2 ) 
cosx pcosm v ' 

sin (/3x) 7i sgn (3 cos (m/3) 
sin x 2 f3 sin m 



o(r 2 ) 



(E.la) 
(E.lb) 



valid as (3 — >■ ±oo with fixed m G (0, tt). 

(IE. lap follows by repeated integrations by parts that bring down inverse powers of 

In (lE.lbp . we split the integral as 



dx 



smx x 



sin(/3x) 



00 , sin(/3x) f°° , sin(/3x) 
dx — — + dx — - 

Jo 



x 



X 



(E.2) 



Repeated integrations by parts now apply to the first two terms in (1E.2j) . and the third 
term equals ^sgn/3 [39]. Combining, we obtain (lE.lbj) . 



F Stationary but non-co-rotating detector 

In this appendix we discuss briefly a detector that is stationary in the exterior region of 
the BTZ black hole but not co-rotating with the horizon. For the transparent boundary 
condition at the infinity, we show that the n = term in the transition rate (14.111) 
breaks the KMS property already in second order in the difference between the horizon 
and detector angular velocities. As the n = term is expected to give the dominant 
contribution when the black hole mass is large, we take this as evidence that the tran- 
sition rate does not satisfy the KMS property, in agreement with the GEMS prediction 

[331 SISSIES]. 
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Consider a detector that is stationary in the exterior region of the BTZ spacetime 
at exterior BTZ coordinate r, but not necessarily co-rotating with the horizon. The 
tangent vector of the trajectory is a linear combination of 8 t and 8$. By ( 14. 3 p and ( I4.4p . 
the lift of the trajectory to AdS3 reads 

Xi = £coshxsinh(2%) , 
T\ = £ cosh x cosh (2ky) , 
X 2 = £sinhxcosh(2|/) , 

T 2 = £sinhxsinh(2?/) , (F.l) 



where we have written sfa = coshx with x > 0, the constant k is proportional to the 
difference of the detector and horizon angular velocities, and y is a parameter along 
the trajectory. We assume \k\ < tanhx, which is the condition for the trajectory to be 
timelike. The proper time r is related by y by r = 2£ sinhxyl — k 2 coth 2 x y. 

Let J rn=0 denote the n = term in the transition rate ( 14. lip . Substituting ( IF.ip 
in (I4.12p . and specialising to the transparent boundary condition, £ = 0, we find 



F n=0 (E) 



1 

4 



1 

2^ 



1 — k 2 coth x 



sin (2E£ sinh v\/l-P coth 2 y y) 
dy V ; A V A J^- . (F.2) 

sinh 2 y — coth 2 ^ smh 2 (ky) 



It can be verified that the quantity under the square root in the denominator is positive 
for < y < oo. 

Expanding ( \F.2\i as a power series in E and then expanding the coefficients as power 
series in k, we find 



1 




4 


+ 




"tt 3 


+ 




12 



7T . 7T ( TX 2 X 
■ — smh Y H 1 , 

4 A 8 V 4 ) sinh X 



cosh ^ 2 + o(^) 



sinh 3 ^ + 



71 



TX' 



sinh x cosh 2 ^ k 2 + O (/c 4 ) 



(^) 3 + 0((^) 5 ) . 

(F.3) 



From (IF. 31) it is seen that the power series expansion of J rn=0 (—E)/J r n= o(E) in E is 
incompatible with a pure exponential in E, and the discrepancy arises in the coefficient 
of the (E£) 3 term in order k 2 . J 7n=0 (|F.2p hence does not satisfy the KMS property at 
small but nonzero k. 
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